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Abstract 

We prove L°°(M. d x [0, oo)) bounds on moments X a := X^meN ma fm( x , t) OI the 
Smoluchowski coagulation equations with diffusion, in any dimension d > 1. If the 
collision propensities a(n, m) of mass n and mass m particles grow more slowly than 
(n + m)(d{n)-\-d(m)) , and the diffusion rate d{-) is non-increasing and satisfies rn~ bl < 
d(m) < m~ b ' 2 for some b\ and 62 satisfying < 62 < b\ < 00, then any weak solution 
satisfies X a G L°°(R d x [0,T]) n L 1 (R d x [0,T]) for every a 6 N and T E (0,oo), 
(provided that certain moments of the initial data are finite). As a consequence, we 
infer that these conditions are sufficient to ensure uniqueness of a weak solution and 
its conservation of mass. 



1 Introduction 

The Smoluchowski coagulation equation is a coupled system of partial differential equations 
that describes the evolving densities of a system of diffusing particles that are prone to 
coagulate in pairs. A sequence of functions f n :R d x [0, 00) — > [0, 00), n G N, is a solution 
of the Smoluchowski coagulation equation if it satisfies 

(1.1) J^»M = d ( n ) A fn(x,t) + Qn(f)(x,t), 
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with Q n (f) = QUf)-Qn(f), where 

n-l 

Qn(f)( x ^) = ^2 a ( m > n - ™>)fm(x,t)fn-m(x,t) 
m=l 

and 

oo 

Qn(f)( x 't) = 2f n (x,t)^2a(n,m)f m (x,t). 

m=l 

We will interpret this solution in weak sense. Namely, we will assume that and Q~ 
belong to L 1 (R d x [0, T]) for each T £ [0, oo) and n G N, and that 

/„(z,t)=S? (n) /°(x) + f s£2Q n {z,8)d8, 

Jo 

where f° denotes the initial data, the semigroup associated with the equation u t = DAu, 
and where Q n (x, s) means Q n (f)(%, s). 

The system (jl.ljl has two sets of parameter values, the sequence d : N — > [0, oo), where 
d(n) denotes the diffusion rate of the Brownian particle of mass n, and the collection a : 
N 2 — > [0, oo), where a(m, n) models the average propensity of particles of masses m and n to 
coagulate. The terms Qn(f) an d Qn(f) are § a i n an< i l° ss terms for the presence of particles 
of mass n that arise from the binary coagulation of particles. 

The system may be augmented by considering the fragmentation of particles into two 
or more sub-particles. A continuous version of the system, in which particles have real (rather 
than integer) mass, has also been considered. The data are defined d : [0, oo) — > [0, oo) and 
a : [0, oo) 2 — > [0, oo), and the sums are replaced integrals in the definitions of Q^if) an d 
Q~(f). The spatially homogeneous version of the equations, in which each /„ is a function 
of time alone, are better understood (the work PQ resolved many of the central questions for 
the discrete case without fragmentation). 

Lang and Nguyen [§] considered a system of mass-bearing diffusing particles, whose 
diffusion rate was chosen to be independent of mass, that are prone to coagulate in pairs 
at close range, and demonstrated that, in a kinetic limit, the density of particles evolves 
macroscopically as a solution of (jl.ip . A kinetic limit means that, with an initial number N 
of particles, the order of the interaction range e = e(N) of any given particle is chosen so that 
a typical particle experiences a rate of collision that is bounded away from zero and infinity 
for all high N. In [7] , jH] , this kinetic limit derivation was extended to permit more general 
diffusion rates, and to include a stochastic mechanism of interaction. Suppose that particles 
of mass iigN have a range of interaction given by r(n)e, for some increasing function r : N — > 
(0, oo). That is, e = e(N) determines the order of the range as a function of the initial particle 
number, whereas the function r specifies the relative interaction range of particles of differing 
masses. The variations to the derivation required by introducing the radial dependence r 
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have been discussed in It is shown there that, if the dimension d > 3, the macroscopic 
collision propensity a : N 2 — > (0, oo) that appears in the Smoluchowski coagulation equation 
satisfied by the macroscopic density profile of the particle system satisfies 

/ \ d-2 

(1.2) a(n,m) < c(d(n) + d(m)) (r(n) + r(m) J , 

where c is a constant without dependence on the details of the stochastic interaction. (Note 
that (3 rather than a was the notation used in [l|,[Hj- The change in notation in this paper 
ensures consistency with the PDE literature.) In fact, the left and right-hand-sides of (jl.2j) 
are of the same order provided that the interaction mechanism is strong enough to ensure 
that a uniformly positive fraction of pairs of particles that come within the interaction range 
coagulate. 

The growth rate of r : N — > (0, oo) presumably depends on the internal structure of the 
particles. If they are simply balls, each of the same density, then r(n) = con 1 ^. On the 
other hand, an internal structure that is fractal might give rise to a relation of the form 
r(n) = n x , for some x — \jd. The physically reasonable range of values of \ would seem 
to be contained in [1 /d, 1] . This is because, whatever fractal structure a particle of mass 
iiGN may have, if contact with another particle is required for the pair to coagulate, then 
the interaction range r(n) is at most of order n, (with the extremal case being that in which 
the particle takes the form of a line segment). Indeed, non-trivial fractal structure would 
suggest that x is strictly less than one. 

Monotonically decreasing choices of the diffusion rates d : N — > [0, oo) seem to be phys- 
ically realistic, since the diffusive motion is presumably stimulated by the bombardment of 
much smaller elements of an ambient gas. The choice d(n) = - in three dimensional space is 
justified if the particles are modelled as balls. (For a mathematical treatment, see j3], which 
derives Brownian motion as the long-term behaviour of a ball being struck by elements in a 
Poisson cloud of point particles of Gaussian velocity.) 

In this paper, we examine the behaviour of solutions of (jl.l)) . We have directed our 
attention to parameter values a and d that seem to be justified by the existing kinetic limit 
derivations of (jl.lj) (though, in principle, other choices may arise from a derivation of the 
equations from a quite different model). 

We mention firstly that, under the assumption 

(1.3) lim — ^ — - = for each neN, 

m^oo m 

the global existence of a weak solution of (jl.lj) has been established in [TO], which work 
includes fragmentation in the equations. From the physically reasonable assumptions that 
d(-) is uniformly bounded and r(n) = o(n), we see that (jl.3j) is satisfied in dimension d = 3 
by choices of a satisfying (jl.2j) . 

An important formal property of solutions of (jl.lj) is mass-conservation. 
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Definition 1.1 Let f solve M.l)) weakly. We say that f conserves mass on the time interval 
[0, T] provided that I(t) = 1(0) for each t G [0, T], where 

I(t) = J^m / f m (x,t)dx. 

While mass conservation holds formally, the only estimate that is readily available is I(t) < 
1(0) for t > 0. Indeed, the inequality may be strict, in which case, gelation is said to occur 
(at the infimum of times at which the inequality is strict). If this happens for a solution of 
(jl.lj) which is obtained as a kinetic limit of a particle system, then, after the gelation time, a 
positive fraction of the mass of particles is contained in particles whose mass is greater than 
some function that grows to infinity as the initial particle number tends to infinity 

In the spatially homogeneous setting, much progress has been made in establishing when 
gelation occurs. In [Q\ and :5j, continuous versions of the spatially homogeneous equations 
with fragmentation are considered. If the continuous analogue a(x,y) of the coagulation 
rates a is supposed to satisfy 

(1.4) a(x, y) = x a y b + x b y a , x, y G R 

with a, b G (0, 1) and a + b > 1, then gelation occurs, unless the inhibiting effect of fragmen- 
tation is strong enough. It is natural to postulate from these results that if the microscopic 
interaction range r(n) we have discussed behaves like r(n) = n x with x > ^Zo, with a diffu- 
sion rate d uniformly bounded below and a bounded domain Q in place of M. d , then gelation 
will occur. This is because the formula (jl.2|) . (which, as already noted, may be written as 
an equality in the case of a reaction mechanism that is not particularly weak), is bounded 
below by the discrete analogue of the coagulation propensity given in (jl.4j) . As we have 
commented, however, we do not anticipate such behaviour for a in equations arising from a 
three-dimensional particle system of the type considered in [TTj . 

Rigorous sufficient conditions for mass-conservation, or for gelation, have been available 
in the spatially inhomogeneous setting only under stringent assumptions on parameters. See 
[3| for the case of constant diffusion rates, and jTS] for a criterion that requires uniform 
boundedness of a and further information about the behaviour of solutions of the system 
(however, each of these papers includes fragmentation in the equations). In Theorem 11.31 
we present a more applicable sufficient condition. 

The result largely depends on new moment bounds on solutions of (jl.lj) . Theorem 11.11 
presents bounds on the L l norm of moments of a solution, and Theorem 11.21 provides L°° 
bounds on such moments. Previously, L°° estimates on solutions of (jl.lj) . with the effect of 
fragmentation included, have been obtained (see [12], [IS], 0,13), under fairly restrictive as- 
sumptions on coefficients for coagulation and fragmentation propensity. The dependence on 
i G N of the bounds obtained on L°° norms of fi does not generally permit deductions about 
the L°°-norm of moments Ylm=i ma fm(x, t) for any a > (note however that such inferences 
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are made in [3] if the diffusion rate is identically constant, or in [12 J if the coagulation rates 
a(n,m) decay quickly enough). 

Our final result, Theorem 1.4, provides a criterion for uniqueness of solutions of (jl.lj) . 
It also relies principally on the moment bounds and is a straightforward adaptation of the 
uniqueness proof of pQ which applies to the homogeneous case. 

Our results for valid for each dimension d > 1. Each deduction, moment bound, mass 
conservation, or uniqueness, depends on some regularity in the initial data, and some as- 
sumption on the parameters of the system. We now state the various assumptions that we 
require. 

Assumption 1.1 

g(n,m) 

ilm 7 w ,/ x u — \T — U - 

n+rrwoo [n + m){d[n) + d(m)) 
More precisely, for every 5 > 0, there exists k = ko(5) > such that if n + m > k , then 

a(n, m) < 5(n + m)(d(n) + d(m)). 

In addition, the function d is uniformly bounded. 

Assumption 1.2 The function d(n) is a non-increasing function of n and that a(n, m) < 
Co(n + m) for a constant Co- Moreover, there exist positive constants r\ and T2 and nonneg- 
ative constants 62 < &i such that, 

r\n~ hl < d(n) < r 2 n~ b ' 2 . 

Assumption 1.3 The function d is uniformly positive and non- increasing, and there 
exists a constant Cq such that 

a(n, m) < Co(n + m). 
Our notation for the various moments of / will be 

(1.5) X a = X a (x,t) = J2^ a fn(x,t), X a = X a (x,t) = Y,™ a d(n) d/2 fn(x,t). 

n n 

and 

Y a (x,t) = ^nm(n a + m a )(d(n) + d(m))f n (x,t)f m (x,t), 

n,m 

(1.6) Y a (x,t) = ^2(n a m + m a n)a(n,m)f n f m . 
We also set 



n.m 



\x\ 2 - d ifd>3, 

;i.7) <j> (x) = { -^log|a;| l(\x\ < 1) if d = 2, 

\x\) l(2|x| < 1) if d= I. 

We now state the four theorems. 
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Theorem 1.1 (Moment bound I) Assume Assumption 1.1. Then for every a > 2 and 
positive A and T, there exists a constant C = C(a,A,T) such that, if 



1.8) J J X a (x, 0)Xi(?/, 0)(f)o(x - y)dxdy < A, 

n,m 



and 

(1.9) ess sup J X a (y, O)0 o (x - y)dy < A, J X a (x, 0)dx < A, 
then 

(1.10) J j ' Y a ^dxdt < C, j J ' Ya^dxdt < C, 
and 

(1.11) sup / X a (x,t)dx < C. 

te[0,T] J 

Moreover, the constant C can be chosen to be independent ofT when d > 2. 

Theorem 1.2 (Moment bound II) Assume Assumption 1.2. Then there exists a function 
7 (a, bi,b 2 ) with lim^oo •y(a,b 1 ,b 2 ) = oo such that if 

(1.12) xyii^iu-^) < °°' x - e Li{Rd x [°' t ])' 

n 

then 

(1.13) £ni/„iu» (R d x[0,T]) < °°> 

/or every e < 7(0, & 1; 6 2 )- 

Remark 1.1 In particular, if Assumptions 1.1—1.2 hold, X a (-,0) E L 1 (M d ) and 

^2n a \\fi\\L°°(»« < oo- 



for every a G N, then X a G L°°(M d x [0, T]) n L^R* 1 x [0, T]) for every a G N and T G (0, 00). 

We refer to (|4.15j) for the explicit form of the function 7. Also note that Theorem 1.1 
offers sufficient conditions to ensure X a E L 1 (R d x [0, T]) (so that this theorem has been 
invoked in Remark 1.1). 
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Theorem 1.3 (Conservation of Mass)Le£ f be a weak solution of (II. Assume that 
()1.10|) holds for a = 2. Then f conserves mass on the time interval [0, T\. Assume instead 
that Xi(-,0) G L°°(R d ), X 2 (-,0) G L 1 (R d ) and Assumption 1.3 holds. Then f conserves 
mass on the time interval [0, oo). 

Theorem 1.4 (Uniqueness) There is a unique weak solution of (jl.lj) on the interval [0, T] 
among those satisfying X 2 G L°°(R d x [0,T]). 

Remark 1.2 Assume that there exist positive constants C\ and c 2 such that ct(n, m) < 
C\[n a + m a ) and din) > c 2 n~ b for all n, m G N. Assume that a + b < 1. As a consequence of 
Theorems 1.1-1.4, if J2n ne \\fn\\L°°(R d ) < °° an d || J2 n ^/nlU 1 ^) < 00 f° r sufficiently large 
e, then (1.1) has a unique solution which is mass conserving. 

Remark 1.3 Theorem 1.1 is also true for the continuous version of the system (jl.lj) with a 
verbatim proof. In the continuous version, all the summations over n and m are replaced with 
integrations with respect to dn and dm. On the other hand, we can derive the continuous 
version of Theorems 1.2 and 1.3 only for a particular solution. The proof of Theorem 1.4 
does not readily adapt to the continuous setting. See Remark 4.1 for further comments. 

If the parameters a and d derived from the kinetic limit of a particle system are such 
that uniqueness among solutions of (jl.lj) is unknown, then, in principle at least, the particle 
system may not approximate a single solution of (jl.lj) . It might, for example, approximate 
several different solutions, each with a positive probability. This would be very peculiar, and 
so, it is pleasing to be able to rule out the possibility by establishing uniqueness. 

If we adopt the relation r(n) = n x for particle interaction range, then, recalling (jl.2j) . 
we have that, in IR d with d > 3, the macroscopic coagulation propensity arising from the 
microscopic random model satisfies 

a(n } m) < c(d{n) + d[mj) max jn, m} x d 

In view of Theorems ll.3l and ll.4j we have a unique solution and mass conservation throughout 
time provided that x £ [0> ^z^) anc ^ suitable moments of the initial densities are finite. As 
such, the discussion following (jl.2j) rules out the occurrence of gelation in three dimensions 
in particle systems of the type considered in [0] and jjj. 

In a similar vein to the comment that follows (jl.4j) . we mention that, in the case where 
d is uniformly positive, Theorem 11.31 gives a sufficient condition for mass conservation of a 
solution of (jl.lj) that is close to being sharp. Indeed, choices of a(n,m) that grow much 
more quickly than Assumption 1.2 permits are bounded below by the expression in (jl.4j) . 
for some choice of a, b G (0, 1) with a + b > 1. Such a choice of a would thus be expected to 
show gelation in the spatially homogeneous case. Corollary 8.2 of i6J adapts the argument of 
the homogeneous case to assert that gelation occurs for any weak solution of the continuous 
version of (jl.lj) in a bounded subset of R d , provided that (jl.4j) is satisfied for such a and b 
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as above. The diffusive motion of particles in M. d may act to inhibit gelation of a solution of 
(II. 1J) . though a dense initial condition is likely to ensure it. 
Acknowledgment. We thank Gabor Pete for comments on a draft version. 



2 The tracer particle approach 

Each of the L°° and moment bounds on solutions of that we present in this paper will 
be proved by PDE methods. However, for each of our results, we earlier derived a similar 
assertion by a quite different approach. A given solution of (jl.lj) is understood in terms of 
the random trajectory of a tracer particle, whose behaviour is typical of the many particles 
that form the density profile of the solution. We have certainly found this random method 
to be intuitively appealing, and it may find application to other PDE for which kinetic limit 
derivations have been made. We have thus devoted this section to explaining the tracer 
particle approach. 

It will in fact be helpful to recall in more detail the model analysed in [7] . A sequence of 
microscopic random models indexed by their initial number iV G N of particles is given. In 
the iV-th model, each of the N particles has an initial location x(0) and integer mass m(0) 
set independently according to 



(2.1) P(m(0) = k) 

with x(0) having law 
(2.2) 



f, 



conditional on m(0) = k. At any given moment of time t G [0, oo), particles of mass k G N 
evolve as Brownian motions with diffusion rate d(k), with d : N — > (0, oo) a given collection 
of constants. Particles are liable to coagulate in pairs when their displacement is of order 
e = e(N) (we refer the reader to the introduction of [7] for the details of the interaction 

mechanism). We set e ~ N" 3 ^ (for d > 3), to ensure that a typical particle experiences 
a rate of collision that is bounded away from zero and oo uniformly in N. At the collision 
event, the two incoming particles disappear, to be replaced by a third, that assumes the sum 
of the masses of the ingoing two, and is located in an e-vicinity of either of the colliding 
particles. While the macroscopic behaviour of the system is likely to be independent of the 
choice of placement of the new particle within this microscopic vicinity of the colliding pair, 
it was convenient for the derivation performed in [7j to assume that, if the masses of the 
colliding pair are n and m, then the location of the new particle is taken to be that of one 
or other of the pair with probabilities —r— and — For what follows, it is convenient to 

r r n+m n+m 
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regard a particle that has a collision as surviving it, and becoming the outgoing particle, 
with a probability proportional to its mass, and disappearing from the model in the other 
event. 

Theorem 1 of [Zj specifies the macroscopic behaviour of this particle system, when the 
initial particle number N is taken to be high. To summarise the result without recourse to 
equations, at typical points (x, t) &M. d x [0, oo) of space-time, the number of particles of given 
mass m G N located in a vicinity of x at time t, normalized appropriately, approximates the 
density f m {x,t), for some solution |/ m : m G N} of (jl.lj) . ^From the result, we may infer 
the law of the trajectory of a typical particle, in the limit of high initial particle number, as 
follows. 

Consider a particle picked uniformly at random from the N particles that are initially 
present. We call this the tracer particle in the iV-th model. In accordance with the preceding 
description, the initial mass m(0) and location x(0) of the particle are given by (|2.1|) and 
Q. We know from [7] that {/„ : R d x [0, oo) [0, oo), n G N} gives the asymptotic particle 
densities of the microscopic models. Hence, asymptotically in high N, the trajectory of the 
particle is such that, if at time t G [0, oo), the particle has location x G M. d and mass m, it 
evolves as a Brownian motion at rate d(m), and experiences collision with a particle of mass 
n at rate a(m,n)f n (x,t). The details of the collision event that we specified imply that the 
tracer particle (x(t),m(t)), on colliding with a mass n particle, survives the collision and 
assumes mass m(t + ) = m{t) + n with probability m ^^ n , and disappears from the model 
otherwise. 

The limit in high N of the law of the tracer particle in the iV-th model is in fact dependent 
only on the given solution {/„ : n G N} that the microscopic particle densities approximate, 
and not on other details of the random models, such as the location of other particles. 
We now formalise this definition: a tracer particle specified in terms of a given solution 
{/ n : n G N} of (jl.lj) and not in terms of the data of any microscopic model. 

Definition 2.1 Let a solution {/„ : n G N} of M-l\) be given. The tracer particle governed 
by f the random process z = (x, m) : [0, oo) — > M d x N U {c} whose initial law is given by 



P(m(0) = m 



with the initial location x(0) having density /„(■)/ J'^d fm{x)dx, conditional on m(0) = m. 
At time t, the particle's location x(t) evolves as a Brownian motion at rate d(m), provided 
that m(t) = m. Moreover, at time t, and for any n G N, the particle is said to undertake 
a mass transition m — > n + m at rate a(n, m)f n (x(t), t) . Such a transition succeeds with 
probability m / (n + m) , in which case, m(t) is set equal to n + m, and fails in the other event, 
in which case, the particle is relegated to the cemetery state c: z(s) is set equal to c for all 
s>t. 
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Remark 2.1 Certain smoothness (or measurability) assumptions are in fact required to 
ensure that the process z = (x, m) : [0, oo) — > M. d x N U {c} exists, even locally in time. 

Let z : [0, oo) -> R d x N U {c} denote the tracer particle governed by a given solution 
{/„ : n G N} of (jl.lj) . Let g n (x,t) denote the density of its location at time t: 



w(x{t) G A,m(t) = n j = J g n (x,t)dx. 

The evolution equation for the system [g n : n G N} is given by 
d 

—g n (x,t) = d(n)Ag n (x,t) 



n-l 



(2.3) + 22 a ( m > n - m )fm(x, t)g n - m (x, t) - 2g n (x, t) a(n, m)f m (x, t). 

m=l m=l 

Note that the choice g n = f n solves this equation. Assuming that there is a unique non-zero 
solution of 1)2.3 J) , we find that 

(2.4) # n = f n . 



Before commenting further on (J2.4)) . we want to emphasise how we have changed our point 
of view of the tracer particle. At first, we regarded it as a typical particle in a microscopic 
random model, whose behaviour in the large is described by a solution of (Jl.lj) . and then, 
secondly, as a random trajectory defined purely in terms of such a solution. The former 
point of view motivates the study of the system (Jl.lj) . We will now discuss how the latter is 
valuable in studying a given solution of (Jl.lj) . Indeed, that the tracer particle defined by a 
given solution of (Jl.lj) may be a useful tool for studying that solution is apparent from 1)2.4)1 : 
if we understand the likely behaviour of the tracer particle, we infer bounds on the density 
{g n : n G N} of its location, and, by ()2.4j) . on the solution {/„ : n G N} itself. 

We have stated the relation (J2.4J) because doing so permits a more succinct summary of 
how the tracer particle approach works. However, in making the approach rigorous, we take 
a different route, which we now summarise. We construct a sequence {f^ : W 1 x [0, oo) — > 
[0, oo),n G N}, indexed by N, which will approximate the density of the tracer particle of 
some solution of 1)1.1)) when iV is high. For each iV G N, the functions for n G N, are 
constructed inductively, on the domains M. d x [0, j^), for each i G N. They are extended to 
R d x [0, ^) by an inductive step in which f^it) for n G N and t G [j^, is defined as the 
density of the location of a tracer particle governed by the constant data f^{i/N) during the 
time interval [-h,t). For each i G N and T G [0, oo), the function is shown to converge 
in L x (lR d x [0, T)) as N — >• oo to a limit f n : R d x [0, oo) — > [0, oo). By mimicking the proof 
of the weak stability result of [10J, we infer that the sequence [f n : n G N} is a solution of 
(jl.lj) . Tracer particle arguments, of which an example will shortly be given, are then applied 
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directly to the tracer particle densities {f^ : n G N}, and the resulting bounds, which hold 
uniformly in N, are inherited in the high N limit by the solution |/ n : n G N}. This method 
for making the tracer particle approach rigorous does have the drawback of applying to only 
one solution of (jl.lj) . 

How may the tracer particle approach be used to prove an L°°-bound on the solution 
{f n : n G N} of (jl.lj) whose construction we have just discussed? To simplify the exposition, 
suppose that the initial condition takes the form J Rd f®(x)dx = 1, f^{x) = for x G M. d and 
m > 1. We will now sketch a proof of the following reformulation of Lemma 14.11 Suppose 
that d : N — > (0, oo) is decreasing, and that d(m) > cm~d( 1 ~ a ) for some a G [0, 1] and c > 0. 
Then 

oo 

(2.5) ^m Q / m (x,t)<C M (x,t), 

m=l 

where u : W 1 x [0, oo) — > [0, oo) solves, 

w = <i(l)A«, 

u(x,0) = / 1 (x,0). 

The solution |/ n : n G N} of (jl.lj) . being given by the density of the tracer particle (j2.4j) . 
has the following interpretation: for iCK^a Borel set, the quantity 

oo „ 

y^rn a I f m (x,t)dx 

is equal to the expected value of the random variable R = R(A) equal to m a if the tracer 
particle at time t has mass m and lies inside the set A. 

In comparing the left- and right-hand-sides of (j2.5j) . we are thus assessing the degree to 
which the dynamics of the tracer particle may increase the expected values of the random 
variables R{A) over those obtained by using a simple Brownian particle (for numerous choices 
of the set A). For example, if the set A is a small ball about x G M d , and the tracer particle 
is close to x at some time s satisfying < s << t, then a mass transition undertaken by the 
tracer particle at times shortly after s will serve to increase the expected value of R, because 
the slower diffusion rate produced by the transition is more likely to leave the particle nearby 
to x at the later time t. However, if a transition occurs that sharply increases the mass of 
the particle, then it is likely to fail, in which case, it contributes zero to the expected value 
of R. This latter effect limits the capacity of the tracer particle to focus towards x. 

Phrasing the question quantitatively, we ask: what is the expected value of R(A) at time 
s < t if a mass transition mi — > m 2 occurs at time s? Assuming that there is no other mass 



11 



transition, and supposing that the tracer particle is at y G M. d at time s, the expected value 
of R is 

1 f f {y-x) 2 



m 



d/2 



(2tt(* - s)d(mi)) 
if no mass transition occurs, whereas, it is 



/n R N m i a 1 /" f (l/-a;) 2 1 , 
(2.D) — m 9 T77T / exp < ; — >dx 

m 2 2 (2n(t-s)d(m 2 )) d/2 Ja ^ 2(t - s)d(m 2 ) J 

if the transition does occur. (The first factor in ()2.6|) is the survival probability for the tran- 
sition). Given that d(m 2 ) < d(mi), the latter exponential term may be bounded pointwise 
by the former, and we find that 



E[ i?if mass transition occurs 



m l d{mx) d ^' 2 m 2 rn\ a d(mi) d ^ 2 



E (R if it does not) m 2 d(m 2 ) d / 2 mf m\ a d(m 2 ) d /* ' 

If a whole sequence of mass transition occurs, rrii — > m i+ i at time U for i e {1, . . . , n}, with 
ti G [0, i] an increasing sequence, we similarly find that the ratio of the expected values of 
R in the case where the sequence of mass transitions occurs and in that where no transition 
takes place, is bounded above by 

1 d(l) d ' 2 
m}~ a d{m) d / 2 

(We have omitted some details of the argument: to make it rigorous, we might use an 
induction on the total number of mass transitions undertaken prior to time t, and invoke the 
strong Markov property.) 

By comparison, the expected value of R if no transition occurs is u(x,t), where 

* = d{1)A% 

u(x, 0) = j\(x, 0). By choosing A = B(x, r) for each r > 0, we learn that 

°° d(l) d/2 
J2m a f m (x,t) < — , u(x,t). 

/,From the hypothesis d{m) > cm~ 2 ^ d ^ l ~ a \ we find that 

°° d(l) d/2 
Y j m a f m {x,t) < -±J— U {x,t), 

m=l 
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as we sought. 

In fact, a more careful tracer particle argument improves this result: the same conclusion 
(12. 5|) may be reached under the weaker assumption that d is decreasing and satisfies d{m) > 
m -(i-a)+e £ or some q, ^ ^ \^ anc [ e g ^ a y This is a reformulation of our claim in Remark 
1.2 which is a consequence of Theorem 1.2. 

We have mentioned that each of our results has an analogue with a derivation that 
considers the tracer particle governed by a given solution of (jl.lj) . We will not outline the 
method of proof of the result on mass conservation (see Theorem 11.3(1 . stating only the 
characterization of gelation in terms of a tracer particle. We alter the definition of the tracer 
particle governed by N}, so that the particle survives every mass transition. Given 

(|2.4J1 . mass conservation of the solution / until a given time T e [0, oo) occurs if and only if 
the tracer particle experiences only finitely many collisions on the time interval [0, T\. 



3 Moment bounds under Assumption 1.1 

Let us first construct an auxiliary function H that will be needed for the proof of Theorem 
11.11 Before doing so, let us make an observation regarding the case d > 3. When d > 3, 
define 



(3.1) H(x) = c{d)\x 



2-d 



where c(d) = (d — 2)~ 1 iu d 1 with ujd denoting the (d — 1) -dimensional measure of the unit 
sphere S ,d_1 = {x G M. d : \x\ = 1}. We then have that AH = —5q, where 5q denotes the 
Dirac's measure at 0. More precisely, for a test function g, the function 

u(x) = J H(x - y)g(y)dy, 

satisfies Am = —g. Note that H > and this property is lacking when d < 2. Because of 
this we can only hope for the existence of a suitable function H such that H > but now 
AH = — So + Error for an Error that can be controlled. This is the content of our first 
lemma. 



Lemma 3.1 Assume d < 2. There exist functions H and K such that H > 0, K is bounded, 
K is of compact support, 

(3.2) -A x H(x) = 5 -K(x), 

and the function H — 0o is bounded. (The function 0o was defined in (jl.7j) . ) 
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Proof. The construction of H for d = 1 is straightforward; we can readily find a nonnegative 
function H such that H = O in [—1/2, 1/2], if = outside [—1,1], and if is smooth off the 
origin. 

For the construction of the function H when d = 2, let us start from the function 0o and 
make an important observation. Note that if 

(3.3) R(x) = J * (f> {x) = — [ \og\x -y\J(y)dy, 

Zn J\x-y\<l 

with J > 0, then R > and 

(3.4) -AR = J-J, 
where 

(3.5) J(x) = — ( J{x-z)dS{z) = J*6 , 

2vr J\ z \ =1 

where dS denotes the Lebesgue measure on the unit circle and Sq denotes the normalized 
Lebesgue measure on the unit circle. In ()3.3|) . we may replace J(y)dy with a measure J(dy). 
Then (J3.4|) is still valid weakly for an obvious interpretation for (|3.5|) . In particular if we 
choose J(dy) = 5 (dy), then 

-A0o = 5 - S . 

Our goal is to replace Sq with a bounded function of compact support. For this we set 
4>i = So * 0o to obtain 4>i> and 

-A(0 O + 0i) = 5o - 5 , 
where Sq = Sq * Sq is now a "function" and is weakly given by 

h(z)So(dz) = /r) _ N2 / / h(z + a)dS(z)dS(a) 

~~ ' |z|=l </[a|=l 
•2tt /-27T 



(2tt) 2 

/ / /i(e iei + j e *)de x de 2 . 

*) Jo Jo 



(27T) 

Note that the Jacobian of the transformation ^ : (# 1; 2 ) l— * (e^ 1 + e* 6 * 2 ) is given by | sin(6 l 1 
6> 2 )|. Since cos(#i — 2 ) — \\e l91 + e t02 \ 2 — 1, we learn that if S (dz) = 5 (z)dz , then 

■SoM = 2 1 L_ =l(|z|<2) 

< 2 ") 2 v'l-dkP-i) 2 

1 2 /t , 

TL(\z\<2). 



(2n)*\ z \y/4-\z\ 2 
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Here the factor 2 comes from the fact that ip maps (exactly) two points to one point because 
^(^lj ^2) = i>(027 The function j is not bounded. Let us apply the above procedure one 
more time to define 4>2 = <Po * So so that 02 > and 

-A(0 O + 0i + ^2) = #o - <?o, 

where 

4 f 1 
(3.6) 5q(z) = 7 —^\ ^ . =l(|z-a| < 2)dS{a). 



(27r) 3 J| a | =1 | z -a|V4 



a 



As we will see, the function 5 is n °t bounded either. However, the function So is less 
"singular" that So- For this we show that in fact So has a logarithmic singularity. To do 
this, first observe that 6o(z) is radially symmetric because both dS and So are rotationally 
invariant. Hence we may assume that z lies on the 37-axis and z > 0. Observe that the 
integrand in (J3.6)) is singular when either \z — a\ = oi \z — a\ = 2. It is not hard to show 
that there exists a positive constant C\ such that if a = e td with 9 G (— 71", 7r], then 

(3.7) ci(|0| + |1 - \z\\) <\z-a\< \6\ + |1 - \z\\. 

The first inequality will be used to treat l/\z — a\ singularity in (|3.6jl . We now turn to the 
singularity which comes from the factor (2 — \z — a|) -1 / 2 . For this observe that if a = e l6> , 
then we have 

7 (0) : = 4 - 1 2 - a| 2 = 3 - z 2 + 2z cos 0. 
Let us choose 60 G [0, 7r] so that 

z 2 - 3 

3 — z 2 + 2z cos 6*0 = or cos 9 = 



2z 

This means that the integrand in 1)3.6)1 is singular at #0 and — 60 ■ Note that when z > 0, the 
singular point #0 exists only if z belongs to the interval [1,3]. Also note that if z is neither 
close to 1 nor 3, then 60 is neither close to nor 7r. From the elementary inequality 

Q _|_ Q Q Q 

7(0) = 7(0) — 7 (#o) — 2z(cos 9 — cos 6*o) = — 4z sin sin , 

we learn that if z is neither close to 1 nor 3, then we can find a positive constant C2 such 
that 



y/A- \z-a\ 2 > c 2 ^\e-e t 







for 6 1 close to 9 . (Recall that z G [1,3].) Since this is an integrable singularity with respect 
to (^-integration, we deduce that So is bounded if z stays away from the circles \z\ = 1 and 
\z\ = 3. We now assume that Izl is close to 1. In this case 9 Q is close to n and \6q — vr I is 



15 



comparable to yfl — \z\\. Also, because of \z — a\ < 2 and 7(6') — 7(6 , ) > we learn that 
9 < 9 . We have 

(3.8) 4 - \z - a\ 2 > c 3 [(tt - 9) + (tt - O )] [(* - 9) - (tt - O )] , 

for a positive constant C3. When |z| is close to 1, the integrand in (jH.fi)) is singular at ±# 
and "almost" singular (see (j3.7)l ) at . From (j3.7)) and (j3.8)l we deduce 

5o(z) < c 4 I log |1 - |z||| , 

whenever \z\ is close to 1. (Here we used the fact that Jq°[(^ — 9) 2 — (tt — 9o) 2 ]~ 1 / 2 d9 is of 
order | log(7r — 9 )\.) 

When \z\ is close to 3, 9 is small (in fact of order 0(y3 — \z\)), and the condition 
\z — a\ < 2 forces 9 G [— 9q, 9q\. On the other hand, 

4 - \z - a\ 2 > c 5 |6> + 9 \\9 - 9 \, 

for a positive constant c 5 , implies that S (z) is bounded for \z\ close to 3. Here we are using 



r9o 

/ (el - 9 2 ) ~ 2 d9 < 00. 

J-9o 



'-Bo 

Putting all the pieces together, we deduce 

(3.9) S {z) < c[|log|l - |z|||] l(\z\ < 3). 
We set 3 = 5 * O to obtain 3 > and 

(3.10) -A(0 O + 0i + 02 + 0s) = S - K 
where 

K(z) = — ( 5 (z + a)dS(a). 

It is straightforward to use (j3.9)l to show that K is uniformly bounded. Now (j3.2j) follows 
from (j3.10j) by choosing if = O + 0i + 02 + 03- It is also straightforward to check that the 
functions 0i, 02, and 03 are bounded. □ 

Let £ be a nonnegative smooth function of compact support with J £ = 1 and set C <5 ( a; ) = 
5~ d ((x/5). We also define = / n * x ( 5 and = * x ( 5 - We certainly have 

(3.11) fi(x,t) = ti(x,0)+ [ t d(n)Af^x,s)ds+ f Q 5 n (x, s)ds. 

Jo Jo 
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Also, as it is well-known, 

(3.12) y^0(n)Q n = ^2a(n,m)((p(n + m) - <f>(n) - (f>{m))f n f m . 

n n,m 

The same identity is valid if we replace f n f m with (f n f m ) * x ( s and Q with Q 5 . Using the 
fact that for <f>(n) = nl(n <£), we have ^2 n (j)(n)Q n < we can readily deduce that 

(3.13) supsupsup / 2, n fn( x , t)dx < oo. 

s e t J , 

n=l 

Lemma 3.2 Lei i7 6e as in Lemma 2.1. Then there exists a constant Co such that 

t t 
sup ^ n Ut *x H) (x, t) < sup Y n Ut *x H) (x, 0) + c t 

x n=l x n=l 

for every positive 5. We may choose Co = when d > 3. 

Proof. We have 

t t „ t t 

^n(fi* x H)(x,t) = Y, n (fn**H)(x,0)- ^nd(n)fl(x,s)ds 

n=l n=l n=l 

,t * /.t « 

+ / y^nd(n)f^* x K(x,s)ds+ l S^nQ 5 n * x H{x, s)ds. 

J ° n=l J ° n=l 

^From the boundedness of K, H > and Y2 n =i n Qn < we deduce 

Y,n(ti* x H)(x,t)<Y,n(fn* x H)(x,0)+ Cl f* f J2nd(n)f*(x, s)dxds. 

71=1 71=1 ^° ^ 71=1 

We now use f!3.13|) to bound the last term to complete the proof. □ 
Proof of Theorem Set 

= f {^n a fi{x^(^nft* x H{x^dx, 

Z ^ = / fXy/n(z,t)) (j^nf n * x H(x,t)j dx. 
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We have that weakly, 

j t Z\t) = - J j (^nd(n)f*(x,t)jdx 

-J fcn a d(n)f«(x,t)\ fcnfi&t)) dx 



.n=l / \n=l 



+ / (j2 na fn^t?j \ J2nQ s n * x H(x,t)J dx 
+ J fen a Q s n (x,t)\ fenf** x H(x,t)\dx 

+ / (jl^fn^t?! (j2 Ud ^ fn*xK( X ,t)^j dx 

+ / (^ nad ( n )fn(x,t)j \^nf°* x K(x,t)j dx 

We now study the various terms which appear on the right-hand side. We certainly have 
Q 1 + Q 2 = -i I y^ l(n,m< £)[n a md(m) + n a d{n)m + m a nd{n) + m a d{m)n]f^f^ n dx 

n,m 

= -- I Hn, m < £)nm{d{n) + d{m)){n a - 1 + m a_1 )/n/£ dx. 

n,m 

^From Yln=i n Qn < 0, we learn that Q 3 < 0. By boundedness of K, (j3.13|) and the bound- 
edness of d(-) we deduce that 

1^5 + ^1 <Ci J X S a dx, 

where X a is defined in ()1.5|) and = X a * ( s . It remains to bound f^. Note that 

i 

^2n a Q n (x,t) = ^[(n + m) a t(n + m<£)-n a l(n<£)-m a l(m<£)}a(n } m)f n f m 

n=l n,m 

(3.14) < c 2 ^2(n a ' 1 m + m a - 1 n)a(n,m)t(n + m<£)f n f m =:c 2 Z 4 . 



^From this, Lemma l3~2"l the boundedness of H — O an d ()1.9|> . 

n A <c 2 J Z{ (j2 n fn ** dx<c 3 J Z*dx, 

where Zf = Z^ * ( s . On the other hand, for 5q > 0, we can find fco = &o(^o) such that if 
&o < then 

Z 4 = ^ l(fco <n + m< C)(n a ~ l m + m a ~ l n)a{n,m)f n f m 

n,m 

(3.15) + m < k )(n a 1 m + m a 1 n)a(n,m)f n f m 

n,m 

< 5 ^ H n + m < £){n a ~ l m + m a-1 n)(n + m)(d(n) + d(m))f n f m 

+ 2k%J2 1 ( n + m < i ) a ( n i m )fnfm- 
n,m 

As a result, 

n,m 

+2c 3 5 [^2l(n,m< i)nm(n a - 1 + m a ~ 1 )(d(n) + d(m))(f n f m ) * ( 5 dx 

n,m 

+2c 3 k% I ^l(n + m < £)a(n,m)(f n f m ) * C^x + c x I X 5 a (x,t)dx. 

n,m 

Here, we are using the identity, valid provided that a > 2, 

(n + m) (n a ~ 2 + m a ~ 2 ) < 2 (n"" 1 + m a4 ) . 
We now send S to to yield 

d -Z{t) < ( 2c 3 5 - ^ y ^ l(n, m < £)nm(n a - 1 + m a - l )(d(n) + d(m))f n f m dx 



dt r w - v J U 2 



n,m 



(3.16) +2c 3 /cq / /J l(n + m < £)a(n, m)f n f m dx + Ci / X a (x,t)dx. 

n,m 

Note that the time integral of the second integral is bounded because 
d ^ 

~J t fn d x < ~ l^2M n + m < £)a(n, m)f n f m dx, 

n=l n,m 
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which implies, 

(3.17) / / '^l(n + m < £)a(n,m)f n f m dxdt < I dx. 

" n.,m n=l 

Furthermore, the equality 

(3.18) - ! J2 n<1 fn(x,t)dx= I J2 na Qn(^t)dx, 



and ()3.14j) imply that 



nt f ^ r ^ rt ns r- 

\ ^2n a f n (x,s)dxds < t / ^ n a f n (x, 0)dx + c 2 / / Z 4 (x,6)dxd6 
Jo J n=1 J n=1 Jo Jo J 

< t [ y2n a f n (x,0)dx + tc 2 [ [ Z 4 (x,6)dxd6. 
J n=1 Jo J 

^From this, ()3.17|) and ()3.16|) we deduce that Z{t) — Z(0) is bounded above by 

(2c3 + C!C 2 t)5 - ^ J J ^2 nmin*- 1 +m a - 1 )(d{n) + d{m))f n f m dxds + c 4 {l + t). 

^ n,m<£ 



We now choose 5o = 5o(t) so that 1/2 > (2c3 + cic 2 t)5o. With this choice, the bounds in 
(fTTUJ) follow. From flTjTUJ) . (l3~T3j) and (I3~TH|) we conclude (fTTTjl . □ 
We end this section with a variant of Theorem 1.1 that holds under Assumption 1.3. 

Lemma 3.3 Under Assumption 1.3, there exists a constant C such that 
(3.19) J X 2 {x,t)dx <(^j X 2 {x,0)dx^j exp ^CT\\ ^ nf°\\ L o 

Proof. We start from 

-g I ^2n 2 f n (x,t)dx = / ^2n 2 Q n (x,t)dx 

n=l n=l 

< / 2 nma(n, m)l(n + m < £)f n fm dx 

n,m 

< f2C j:n m{ n + m )l(n + m < e )f n f mix 

n,m 



dx. 
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This and Grownwall's inequality imply (|3.19J) because we can use the uniform positivity of 
d(-) and Lemma 3.1 of Section 3 to assert that X\ G L°°. □ 



4 Moment bounds when d(-) is non-increasing 

This section is devoted to the proof of Theorem 11.21 We start with a lemma. 
Lemma 4.1 Assume d(-) is non-increasing. Then 

oo 

(4.1) = J2 nd ( n ) d/2 fn(x,t) < d{l) d / 2 U {x,t). 

n=l 

where u is the unique solution to ut = d(l)Au subject to the initial condition u(x, 0) = 
J2n=i n fn(x,0). 

Proof. We first establish 

(4.2) j^nd{n)^f n {t) < d(lf 2 S^ (l>/s) + d(£) d / 2 JjS (f^nQ^ ds. 

Here for simplicity, we do not display the dependence on the x-variable. Note that (|4.2|) 
implies 



(4.3) ^{n)^f n < d(l)^S? W ( n fn 



because Yli n Qn < 0. Evidently ()4.3|) implies (|4.1j) . 

We establish (14. 2|) by induction. ()4.2|) is obvious when £ = 1 by definition; in fact we 
have equality. Suppose (|4.2J) is valid. We would like to deduce (I4.2J1 with £ replaced with 
£ + 1. To do so, first observe that if D\ > /} 2 and g> > 0, then 

(4.4) D/^S^g > D 2 d/2 S? 2 g. 
^From this and (|4.2jl we learn 

(4.5) ^nrfH rf / 2 /n<rf(l)^ (1) ("tnA + d(£ + if 2 J* S^ s +1) fcnQ n (s)\ ds 
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because d(£) > d(£ + 1) and Y?i n Qn < 0- Applying (JO) to 

f i+1 (t) = Sf +1) f l+1 + f StT ] Qt + i(s)ds, 

Jo 

yields 

(4.6) u +l {t) < (^^)" /2 ^V,% + ^S +1) g, + i( s )^. 



We multiply both sides of by (£+ l) d/2 and add the result to (jPj) . The outcome 

is 



This completes the proof. □ 

Remark 4.1. For the continuous version of (jl.lj) . a similar proof leads to an L°° bound 
on J °° nd{n) d l 2 f n dn provided that the function d(-) is piecewise constant, uniformly positive 
and nonincreasing. If we assume only the second and third conditions on d(-), we can 
establish the same bound on a solution provided that this solution can be approximated 
by solutions corresponding to piecewise constant d(-). Of course, if we already know the 
uniqueness of solutions to (1-1), then our L°° bound applies to all solutions. But we do not: 
when uniqueness is proved in Section 4, we will use Lemma 3.1 and its consequence Theorem 
1.2. 

Proof of Theorem \1.2L 

Step 1. Let us simply write L l for L e (W d x [0,T]). We first show that if I > 1, then 

(4.7) X ae+1 G L 1 X a+1 G L l . 

Indeed, 

nd(n) d / 2 f n 



n 



by Holder's inequality. This implies ()4.7|) because by Lemma 3.1, X\ G L c 
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Recall that we are assuming 

r\rr bl < d(n) < r 2 n~ b2 . 
^From this and (|4.7|l we deduce that 



This means that 

2a+b 2 d — 2 

(4.8) x a e L 1 => x b e L 2 w- 2 , 

provided that 2b + b x d - 2 > and (2a + b 2 d - 2) /(26 + - 2) > 1. 

Step 2. We now try to bound with the aid of (|4.8j) . Using a(n, m) < C (n + m), we 
certainly have 

Qn = Yl a ( n l' n 2)/m/n 2 

n 1 +n 2 =n 

< ^2 M n i ^ n / 2 or n 2 > n/2)a(ni,n 2 )fmfi 
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ni,ri2 

< 2C [X 1 X (n/2) + X l (n/2)X(0)} 
where X a (iV) = Y, m >N ma fm- Hence 

< c 2 n- £ [X 1 X e + X 1+e X ], 

WQ+Wl* < ^[IIXxII^iix.ii^ + iix^ii^IIXoII^] 

provided that \ = ji + j^ = j^ + ji- To use (|4.8j) . let us first assume that b\d> 2 and that 
2a + — 2 > &id + 21. Assume X a e L 1 . Then we use (14.8(1 to assert that 

2a+b 2 d-2 2a+b 2 d-2 

X 1 e L , X G L <>i d - 2 

Hence if we set 

2a + b 2 d-2 



(49) P 2 0l d + 2£-2 : 

and assume that p > 1, &id > 2, then we have that 

(4.10) X a 6 L 1 => HQji, < 



c n 
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However, if bid < 2, then we have that Xq £ L°° because Ai e L°° . ^From this we learn 
that in this case (|4.9jl is true but now for 

(4.11) 2a + M-2 

1 J P 2t + b x d 

Step 3. Note that if 



p D {x,t) 



(47r^)- d / 2 exp(-g) iff>0, 
if t < 0, 



and g is a function with g(x,t) = for t < 0, then J S^_ s g(x, s)ds = (pr> * g)(x,t) where 
the convolution is in both x and t variables. Also note that 



T r rT r / i \dr , ^ 





{po{x,t)) r dxdt = I I ( —j=\ pi (—==,! ) dxdt 



= f {Dt)i^dt = c(T, r)D^-^ 
Jo 

with c(T, r) < oo if and only if r < | + 1. We certainly have 

/„(*,*) = (^ (n) /„°)(^) + ( Pd (n) * Q„)(X,0- 

So, 

(4-12) ||/„|Uoo < \\(S? n) fZ)(x)\\ L - + \\p d(n) \\ L r\\Q n \\ LP 

provided that i + - = 1. Since pz> £ L r with r < | + 1, it suffices to have 

1 1 

+ - < 1. 



l + 2/d p 
Choosing p as in (|4.9|) or ()4.11|) requires 



2bid + 2£-2 b x d + 2£ 2 
or — < 



2a + b 2 d-2 2a + b 2 d-2 d + 2' 
More precisely, 

(-^(2a + b 2 d-2)-b 1 d+l if M > 2, 

' [^(2a + M- 2) - ±M ifM<2. 
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In summary, we need £ to satisfy ()4.13|) and r to satisfy ~(1 — r)>— 1. As a result, if 
X a G L 1 and £ satisfies (|4.13jh then 



where 



H/nlU- <A n + c 3 n- £ d(n) d ^/ 2 <A n + c 4 n"^)" 1 

^=ll(^ (n) )/n (^)ll^<||/„ |Uoc (Kd ). 



Final Step. We have that ^ n ^ e ||/nlU°° < oo if X a G L 1 , 



(4.14) 
and, 

For (|4.14jl it suffices to have 

In other words, 

(4.15) e< 7(0,61,63) := 



< 00, 



< OO. 



6l < -1. 



^(2a + 6 2 d - 2) - 6i(d + 1) ifM>2, 



^(2a + 6 2 (i-2) - |6i(i-6i -1 ifM<2. 



□ 



5 Uniqueness 

The main result of this section is Theorem 15.11 Theorem 1.4 is an immediate consequence 
of Theorem 15.11 

Theorem 5.1 Assume that a(n,m) < c^nm and let f and g be two solutions with 



where W abbreviates L p {R d x [0,T]). Then 



n2 9r. 



< A, 



x (t) '■= J ^2 n \fn - g n \(x,t)dx 
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satisfies 

(5.1) X(t) < e 4coA *X(0), 

fort < T. In particular, «// n (-,0) = g n (-,0) for all n, then f n {-,t) = g n (-,t) for all n and 

te[o,T). 

We first state a straightforward lemma: 
Lemma 5.1 Let u be a weak solution of 

u t = DAu + h 

with u and h e L l . Assume that ip is a continuously differentiable convex function with 
\ijj'(a)\ < ci for a constant c\ and all aGl. Then 

(5.2) J ip(u(x,t))dx < J ip(u(x, s))dx + J J ip' (u(x,6))h(x,6)dxd9 
whenever < s < t. 

Lemma 15. II is established by choosing a smooth mollifier p e , and showing the inequality 
(5.2) for u t = u* p e and h e = h* p e . We then pass to the limit e — > 0. We omit the details. 

Proof of Theorem 1^.1 . Choose a continuously differentiable convex function ips so that 
ip' s (r) = sgn(r) for r (-5,5), ips e C 2 , \ip' s (r)\ < 1, ^(O) = and ips > 0. We then apply 
Lemma f5. II to assert that the expression 



(5.3) / Y, n Mfn(x,t) - g n (x,t))dx, 

J n=l 

is bounded above by 

N 

^ n if) S {fn(x, s) - g„(x, s))ds 

n=l 

+ / / Z>^(/nM) -g n (x,0))(Q n (f)(x,0) -Q n (g)(x,6))dxd6 

Js J n=l 

N 

= / /^ipsjfnjx, s) - g n (x, s))dx 

J n=l 

m Tn ^m)(fnfm gngm)dxd9. 



N 

(5.4) 
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where T n = n ip' s (f n — 9n)^(n < N). Observe that if n, m < N, then 

(r«+m Tn ^m)(,fnfm 9n9m) i 

equals 

(J^n+m -T n ^m)\fn 9n) fm (J^n+m Tji ^m)\fm 9m) 9n 
< (n + m)\f n - 9n\fm ~ nip' 5 (f n - 9n){fn ~ 9n)fm + m\fn - 9n\fm 

+ (n + m)\f m - g m \gn + n\f m - g m \g n - mijj' s (f m - g m )(f m - g m )g n 
< 2m\f n - g n \f m + 2nl(\f n - g n \ < 5)\f n - g n \f m 
+2n\f m - g m \g n + 2ml(\f m - g rn \ < 5)\f m - g m \g n . 

/From this, (5.4) and a{n,m) < conm, we learn that the expression (5.3) is bounded above 
by 

N 



5^ ni) & {f n {x, s) - g n (x, s))dx 



n=l 



+2c Vl 



+2c S 




N 



^2n\f n -g„ 




n=l 
N 



N 



^2m 2 (f m + g m ) 



,m=l 



dxdd 



J2n 2 M\fn-9n\<5) 



n=l 



N 



^2 + 9m) 



rn=l 



dxdO. 



We then send <5 ^ and N — > oo, in this order to obtain 

/oo „ oo 

^n\f n (x,t) - g n (x,t)\dx < / ^ n\f n (x, s) - g n (x, s)\dx 
n=l J n=l 



+2c 




n=l 



9m) 



m=l 



dxdO. 



The theorem now follows from this and Gronwall's inequality. 



□ 
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6 Mass conservation 

Proof of Theorem IJ..VL We first assume that Y\ G L 1 . Evidently, 

J f n f n \ dx = — < N < n + m)n + l(m < N < n + m)m} 

\n=l / n,m 

= —2 l(n < A" < n + m)na(n, m)f n f m 

n,m 

> -2 l(n > iV/2 or m > N/2)nma(n, m)f n f m 



d_ 

dt 



a(n,m)f n f n 



The limit iV — > oo of the time average of the right-hand side is because Y\ G L 1 . From 
this we can readily deduce that 



lim 



j ( ^2 n fn{ x ^)J dx ~ J ( ^2 n fn{x, 0) j dx 



0. 



This completes the proof when Y\ G L 1 . 

We now assume that Assumption 1.3 holds. We have, 



d 
dl 



/ 1 \ dx = — < A" < n + m)n + l(m < N < n + m)m}a(n, m)f n f„ 

\n=l / n,m 

> -2C ^2l(n<N<n + m)n(n + m)f n f m 

n,m 

> -2C M n < N/2, m > N/2)n{n + m)f n f m 

n,m 

-2C H n > N/2)n(n + m)f n f m 

— — — 
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We certainly have, 





C 










n l 1 


m>N/2 


+ 


c 










n L 2 


m>7V/2 


^2 || X 1 < 


c 




E m ^ m 






n>N/2 


j2 m 


+ 


c 




^ ^ /m 






n>N/2 


Ll m L 



where L p abbreviates L p {R d x [0,T]). Since Y. m m fm e £ 2 , E m m Vm G L\ and En n /« G 
by Lemmas 4.1 and 2.3, we are done. □ 
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